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Kształt i rozmiary jąder o 50 < Z, N « 80 


INTRODUCTION 


The neutron deficient nuclei with 50 « Z, N « 80 nucleon numbers 
have been already investigated theoretically [1,2] but the last development 
of the laser measurement techniques [3,4] as well as the new possibilities of 
the microscopic analysis enable now the better explanation of the diversing 
sizes and shapes mechanisms in this region. The Te-Gd isotopes far from 
the 8 stability line show some interesting features. They are well deformed, 
with the deformation energy up to ~ 10 MeV. Their potential energy reaches 
minimal values not only for the prolate shapes but also for the oblate ones, 
especially for isotopes with N > 74. They have also some hexadecapole 
deformation at the equilibrium point [5]. 

For a long time there was a shortage of the single-particle levels scheme 
parameters here. The extrapolated from the rare earth and actinide regions, 
Nilsson potential parameters, were not sufficient to describe the subtle effects 
in the neutron deficient nuclei. Now we have used the universal set of 
parameters [6] depending only on the average mass number of the whole 
region: A ~ 126 in our case. l 

We have performed the dynamical calculation on the basis of the 
collective hamiltonian, obtained in the generator coordinate method (GCM). 
This hamiltonian consists of the GCM mass parameters and the potential 
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energy improved by the zero point correction term. The many-body 
hamiltonian used here contains in spite of the Nilsson type single particle 
hamiltonian the pairing interaction and the long range two-body forces in 
the local approximation [7]. The quadrupole (e) and hexadecapole (€4) 
deformation parameters are taken as generator coordinates, and BCS ground 
state as the generator function. The BCS wave function was approximately 
projected on a good particle number [8] and the average strength of pairing 
forces was taken from ref. [9]. 

We have calculated the potential energy surfaces, equilibrium deforma- 
tions, mass parameters and dynamical multipole moments (A = 0,2,4) 
for all the even-even nuclei with nucleon numbers 50 « Z « 80, Z € N « 80. 
The chief points of the theory are presented in chapter 1. The numerical re- 
sults of the calculation are illustrated in section 2. As the isotopic shifts of 
the mean square radius do not agree with the experimental data the more 
accurate investigation of their pairing and deformation parameters depen- 
dence is presented in chapter 3. The conclusions and proposals of further 
investigations are drawn in the last part of the paper. 


THEORETICAL MODEL 


The calculation was done on the two dimensional grid of deformation 
parameters —0.4 € ez = € € 0.5, —0.12 < «4 < 0.12. 

The single particle (hamiltonian H. p) eigen problem was solved with the 
Nilsson potential with the new correction term parameters proposed by Seo 
in [6] ko = 0.021, xı = 0.90, yo = 0.062 


Hspln) = enin). (1) 


_ The many body hamiltonian H consists of the mean field hamiltonian 
Ho of Nilsson type taken in the given grid points fe,) = (e,e4), the pairing 
forces and the long range two body correlations in local approximation [10] 


B+ olfa}- Y; xf UNT, — (2) 
Hw=2,4 
where 


Ho = Y eC C, (3) 
n 


with the Cy, Cs fermion annihilation and creation operators. 
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The strength of the long range forces xy is obtained from the selfconsi- 
stency condition 


a_ O(a| F. | a) 


Xu "TRB > Hea) - (4) 
The F, operators are taken in the form 
oH oH 
Elta) = la -lel ga 14) la - (5) 


The operator $ is 1 
SERGE 
n 


The pairing strengths for protons (Gy) and neutrons (Gn) are equal 


GZ?! = G,N?/5 = 0.29 Wo . (6) 


The eigen-function of the hamiltonian (2) is approximated by the BCS wave 
function depending on the single particle coordinates (z) and parametrically 
on the collective ones a: 


la) = la, (23) = Į [(Un VC), (7) 


where V2 is the pair occupation probability, UZ = 1 V/, |0) is the particle 
vacuum state. The function |a) will be used as a generator function and the 
deformation parameters are taken as the collective generator coordinates 
Qj, = €,02 = €4. 

The collective hamiltonian Hon is given in the generator coordinate 
method as 


Recon = 7 +0. (8) 
The kinetic term 7 is 


T=- b TET ja; S GB) 5 De (9) 


where the mass parameters are 
gy zy ha) (10) 


The tensor overlap width is 


— = 


00 
Jij = (a| da; 0a; la) (11) 
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and jh, is a linked matrix element of the many-body hamiltonian 


— 


o 4 7 is 
jhk = (a da; H rim la) — (a| H Ja) . (12) 


The potential term of the collective hamiltonian V consists of the two parts 
Y = (alla) — Eo (13) 
where Eo, so called zero point correction is equal here 
E i os 
Eo = * Bh; (14) 


The expectation value of the many body effective hamiltonian (2) is 
evaluated within the Strutinsky [11] prescription ( EsrRuT) 


(a ja) & Estrur (15) 


The Strutinsky energy reproduces well the nuclear masses and consists of the 
macroscopic liquid droplet ELp [12] part and the shell correction AESHELI. 
describing the shell and pairing effects on the potential energy 


Estrur = ELD = AESHELL (16) 


The eigen problem of the collective hamiltonian Heott (7) is solved by 
diagonalisation in the two dimensional harmonic oscillator wave functions 
base 


— 


Heon®a * tata (17) 


The full many-body wave function Ya describing a nucleus in a state a is 
given within the GCM approximation by the integral: 


val{z}) = | fala)la, {2})da (18) 


This function will serve to calculate the dynamical values of the multipole 
moments. The weight function fa (a) in (18) is directly connected [13] with 
the collective wave function 9, from eq. (17). 

The multipole moments operators are defined as follows 


Qo = 1?,Q2 = r?Pa(cos 9), Q4 = r*Pa(cos b), (19) 


where r, V are the single particle coordinates, P, — the Legendre polyno- 
mials. i 
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It was proved in [13] that the expectation value of Q, operator between 
the ground state functions o is equal to 


Qa = | öülalõnla)õoda, (20) 


where the integral is evaluated in the collective space only. 

In order to compare the theoretical estimates of the quadrupole moments 
to the experimental data, obtained from the reduced quadrupole transition 
probabilities B(E2,2* — 0*) [14] should use the sum rule 


' 
5 a 
+ — ot) = — 2 
RA 0*)= 167 f BoQ?$oda (21) 


where the sum goes over all the possible 2* states and is model independent. 
There is usually one transition B2 favoured in this sum [15] so it is reasonable 
to compare the available experimental data obtained from B( £2) transitions 
[14] with the quadrupole moments calculated as follows 


(93^ = [f eee, (22) 


The mean square radii of nuclei (r?) are usually not given in experiment 
straighty, we deal rather with the isotopic shifts of (r?) between various 
mass numbers [3,4] 


Ar) = (I (23) 
They are related to the electric monopole moments Qo by 


Foe e Qo/Z (24) 


If the mean square radius for the magic spherical nucleus was unmeasured 
we took its liquid drop estimate, designing the experimental OO up to this 
constant. 


2. THE EQUILIBRIUM DEFORMATIONS, 
POTENTIAL ENERGIES AND MOMENTS 


In the table 1 are listed: the equilibrium values of the deformation 
parameters e, eg, potential energies V and deformation energies 


Ede = ve, e4) m V(0, 0), (25) 
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Then the electric monopole moments Qo and the isotopic shifts of mean 
square radius 
alr?) - (DATE P GE (26) 


are compared with their available [3,4] experimental values 0(r?)*. The 
electric quadrupole moments are calculated as (Q2) and (Q 02)? The second 
ones are closer to experimental data [14] . The corresponding reduced 
quadrupole transitions 


B2 = iQ» (27) 


52 R is: (02 (28) 


and the experimental values B2**? are also presented. The hexadecapole 
moments Q4 can be seen in the last column of the Table 1. 

The results are printed for all the even-even combinations of 50 « Z « 80 
and Z « N « 80. 

For the most isotopes the prolate equilibrium shapes are favoured. Only 
the Te nuclei and the heaviest isotopes of other elements have a chance 
to be oblate. The deformation energies, thanks the influence of zero point 
vibration to became larger by up to 1.5 MeV in our calculation, and the far 
from magic numbers nuclei are usually very well deformed. It means that 
the nonaxial y instability in this region would be not so important as it was 
suggested in ref. [1]. 

Presented in Fig. 1 the mean square radius (points) calculated in fact 
up to the value of the spherical lightest isotopes radius do not reproduce 
well the data (crosses) in our calculation, specially their isotopic shifts show 
sometimes quite opposite behaviour to the experimental data. For example 
the theoretical shifts of the heaviest Te, Xe, Da isotopes grow with A while 
their experimental values decrease. The more detailed analysis of our model 
is needed and will be presented in the next chapter. 

In Fig. 2 [16] the quadrupole moments are drawn. The agreement 
with experimental data (crosses) of the (Õ2)1/? values (circules) is in 
general better than (Oz) (points) specially for Te and heavy Xe, Ba, Ce 
isotopes. The negative moments (O2) are obtained for the heaviest Ce-Gd 
isotopes. 

The hexadecapole moments are shown in Fig. 3. The systematic smaller 
values are obtained for N = 70, because of the shell effects. 

In general the shapes of the nuclei are well described but their sizes 
demand the more careful description. 
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66 78 58 66 . 78 62 78 70 8% A 


Pre PFT TT TITAN 


?)(b) 
» (rzy T, 
* exp 


Te Xe Ba Ce Nd SmGd 


LALA 1 LL: LILLA 


118 130 112 120 132 118 


Fig. 1. The microscopic electric mean square radii (r?) in b (points) compared with 
the experimental data [4] (crosses) for even-even Te-Gd isotopes 


78 62 78 70 E. * * 
eee TT 
Q,lb) 

Te 


-4 
E obit ttl 
106 108 130 112 120 132 118 134 126 136 30 38 136 BB 
Fig. 2. The electric quadrupole moments Qz in b calculated as (02) (points) and 
(Q3) (circles) compared with experimental data [14] (crosses) for even-even Te-Gd 
isotopes 
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Fig. 3. The microscopic electric hexadecapole moments Q, in b? for even-even Te-Gd 
isotopes 


3. THE DEPENDENCE OF MEAN SQUARE CHARGE RADII 
ON THE MODEL PARAMETERS 


In order to improve the reproduction of the mean square radius we 
analyse their dependence on the diversing parameters [17]. 

The presented in Fig. 1 experimental data [4] are obtained assuming 
the radius of the magic isotope !3?Ba and !4°Xe from spherical liquid drop 
model. It is seen that the microscopic (r?) grow quicker with A then 
the experimental ones, even when their average values are close to each 
other. 

'The mean square radii depend mostly on the quadrupole deformation e. 
In Figs. 4 a,b we can see this dependence around equilibrium deformation 
€) marked by the arrows for Xe and Ba isotopes. Also the hexadecapole 
. deformation influences the (r?) much as one can see on Figs. 5 a,b for Xe, 
Ba isotopes. Note the almost one order of magnitude difference in the 
(r?) scale of Figs. 4 and 5, 6. Because these collective variables are not 
sufficient to give a good behaviour of 6(r?) we should take next the pairing 
— gap parameter A as a dynamical (collective) parameter. As it is shown 
in Figs. 6 a,b the (r?) value of stable Ke or 12PBa nuclei grows strongly 
with A. 

It is almost impossible to get the values of (r?) and Q} — close to the 
experimental data including the e, «4 parameters only. We have evaluated 
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Fig. 4. The dependence of the microscopic (r?) in fm? on the quadrupole deformation e 
for the Xe (a) and Ba (b) nuclei. The equilibrium deformations e” are signed by 
arrows 
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Fig. 5. The dependence of the microscopic (r?) in fm? on the hexadecapole deformation 
e for the Xe (a) and Ba (b) nuclei. The equilibrium deformations e? are signed by 
arrows 
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Fig. 6. The dependence of the microscopic (r?) in fm? on the pairing energy gap A (in 
MeV) for the Xe (a) and Ba (b) 


within the liquid drop model the values of (€' PAP) (squares in Fig. 7 a,b) 
which would reproduce the experimental data 


(6, rad 9 (45, 4) £ (0,8) 


We have assumed the favored by the microscopic calculation prolate 
deformation of the nucleus when solving the above equation. The solutions 
are far from the microscopic (eb, eg) (stars in Figs. 7 a, b) equilibrium 
deformations both for Xe and Ba isotopes, as it is seen in Figs. 7 a,b 
respectively. 

Let us try to reproduce the experimental data (r*)*P,Q3'P with the 
quadrupole deformation e and the energy gap A varying them free within 
the microscopic BCS model 


((c, Or) BSS (LD A) z (e, A9) £ (e, A") 


As you can see in Figs. 8 a,b the extracted from experimental ((r?) e, 
Or), (erb, A) values (squares) do not agree either with the microscopic 
(e, A?) (stars) equilibrium points, nor even with the values (e, A”) 
obtained from the experimental masses [18] of nuclei (circles). The values 
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Fig. 7. The microscopic equilibrium points (ee, ef) (stars) and their estimates deduced 
from the experimental values of ((r?) *, OY) [3,4,14] within the liquid drop model 
(erb, LD) (squares) for the Xe (a) and Ba (b) isotopes 


of A which reproduce the experimental (r?)*7P and QS” are almost half 
of their estimates from masses A" It is necessary to use the dynamical 
pairing model [19] which would give the most probable A values smaller 
than the equilibrium value corresponding to the minimal potential energy. 
Such calculation is prepared now. 


CONCLUSION 


The following conclusions can be drawn from our calculations: 

1. The dynamical treatment should be used in order to describe the 
nuclear shapes. 

2. The quadrupole moments should be calculated as (63)? rather than 
as (02) in order to compare them with the data experimental transition 
probabilities. 

3. The nuclear radius description for the far from the £ stability line 
nuclei demands the dynamical treatment of pairing forces, especially the 
inclusion of the energy gap as a dynamical variable. 
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Fig. 8. The microscopic values of (e, A?) (stars) compared with the evaluated from exper- 
imental masses (e, A”) (circles) and the (eP, A) (squares) reproducing in microscopic 


model the experimental ((r?)**P, O) for Xe (a) and Ba (b) isotopes 
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STRESZCZENIE 


Przeprowadzono mikroskopowe dynamiczne rachunki elektrycznych momentów o po- 
lowosci A = 0,2,4 parzysto-parzystych jąder neutrono-deficytowych. Zaproponowano 
nową interpretację teoretyczną doświadczalnych momentów kwadrupolowych otrzymywa- 
nych ze zredukowanych prawdopodobieństw przejść B(E2,2* -> 0*). Przeprowadzono 
dokładniejsze badania zależności kwadratów średnich promieni jąder od parametrów de- 
formacji wyższych multipolowości i oddziaływania „pairing”. 
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